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ABSTRACT 

Highly precise and nearly uninterrupted optical photometry of the RR Lyrae star AQ 
Leo was obtained with the MOST (Microvariability & Oscillations of STars) satellite 
over 34.4 days in February-March 2005. AQ Leo was the first known double-mode 
RR Lyrae pulsator (RRd star). Three decades after its discovery, MOST observa- 
tions have revealed that AQ Leo oscillates with at least 42 frequencies, of which 32 
are linear combinations (up to the sixth order) of the radial fundamental mode and 
its first overtone. Evidence for period changes of these modes is found in the data. 
The other intrinsic frequencies may represent an additional nonradial pulsation mode 
and its harmonics (plus linear combinations) which warrant theoretical modeling. The 
unprecedented number of frequencies detected with amplitudes down to millimag pre- 
cision also presents an opportunity to test nonlinear theories of mode growth and 
saturation in RR Lyrae pulsators. 

Key words: stars: variables: other - stars: pulsating - stars: individual: AQ Leo - 
photometry: spacebased - methods: data analysis - methods: statistical 



1 INTRODUCTION 

The human part of the history of AQ Leo (V ~ 12.6) be- 
gins with its classification as a variable of the RW Auri- 
gae type by Hoffmeister (1944) at the Sonneberg Observa- 
tory. Fiftteen years later, Wenzel (1961) concluded that the 
star is an RR Lyrae pulsator with a strongly variable light 
curve. Jerzykiewicz & Wenzel (1977) thoroughly analysed 
their photometry gathered during 1973 - 1975 to recognise 
for the first time that AQ Leo is a double-mode RR Lyrae 
star. They identified the prototype of the class. Jerzykiewicz, 
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Schult & Wenzel (1982) later re-evaluated the data and their 
comparisons to earlier archived exposures raised the possi- 
bility of a period change in the early 1970s. In the meantime, 
double-mode RR Lyraes, designated RRd stars, have been 
identified in globular clusters, dwarf spheroidal galaxies, the 
Magellanic Clouds and the Galactic field (see, e.g., Clement 
et al. 1991). 

RRd stars pulsate in a combination of first overtone 
and fundamental radial modes. The ratio of the first over- 
tone to the fundamental period is about 0.74 - 0.75 (see, 
e.g., Alcock et al. 2000). Usually the amplitude of the first 
overtone is stronger, but not always so (e.g., Oaster et al. 
2006). RRd variables are particularly important because 
their double-mode nature affords an opportunity to deter- 
mine their masses based on the period ratio, largely inde- 
pendently of stellar evolution theory, and hence to a much 
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higher precision than is possible for monomode RR Lyrae 
stars (Kovacs et al 1991). 

AQ Leo happens to fall in the field around a MOST 
space mission Primary Science Target, i Leo. This align- 
ment presented an ideal opportunity to obtain the first new 
light curve of this star in three decades, with very com- 
plete time coverage and unequalled photometric precision. 
The long time baseline between the JW77 photometry and 
the MOST observations is well suited to search for period 
changes in AQ Leo. The precision of the MOST photome- 
try means that the amplitude ratio of the two modes can 
be measured more accurately than before; however, in the 
MOST custom passband, not in a standard filter system. 
If AQ Leo is reobserved by MOST in the future, it will be 
possible to test theoretical predictions of mode switching in 
RRd stars. The sampling and precision of the MOST pho- 
tometry also make it possible to search for nonradial modes 
in an RRd star with a sensitivity that cannot be achieved 
in groundbased data. 



2 MOST PHOTOMETRY 

The MOST (Microvariability & Oscillations of STars) space 
mission was designed to obtain high-precision optical pho- 
tometry of bright stars, monitored for weeks at a time with 
high duty cycle and time sampling of at least once per 
minute (Walker et al. 2003). MOST is a microsatellite hous- 
ing a 15-cm Maksutov telescope feeding a CCD photometer 
through a single custom broadband filter (350 - 700 nm). 
Its polar Sun-synchronous orbit allows it to observe stars in 
its Continuous Viewing Zone for up to two months without 
interruption. 

AQ Leo was observed by MOST satellite in the same 
field as its Primary Science Target t Leo, for 34.434 days 
during 14 February - 22 March 2005 (HJD 2453416.671 - 
2453451.105). One of the faintest stars ever observed by 
MOST at V ~ 12.6, AQ Leo was a Direct Imaging tar- 
get. Direct Imaging is the MOST observing mode where a 
defocused star image is projected onto an open area of the 
Science CCD (see Rowe et al. 2006). The star is sampled 
within a 20 x 20-pixel subraster of the CCD and photome- 
try is extracted through aperture methods and PSF (Point 
Spread Function) fitting of the image profile. The exposure 
time was 25 s and exposures were made every 30 s. The data 
were reduced independently by two of the authors, JFR at 
UBC and DH at IfA Vienna, and the results then compared. 

The UBC MOST Direct Imaging photometry reduction 
pipeline is similar to the treatment of groundbased CCD 
photometry in that it measures the stellar flux though a 
combination of an aperture of pixels centred on the star and 
fitting of its Point Spread Function (PSF) image profile. To 
minimise the effects of pointing errors on the photometry, 
dark and flatfield corrections were performed by monitoring 
individual pixel responses during test exposures on fields 
empty of stars brighter than the background. The correla- 
tion in the raw photometry between the instrumental mag- 
nitude light curve and the estimated sky background was 
removed as described in Rowe et al. (2006). 

The IfA Vienna reduction was based on a pixel-to-pixel 
decorrelation technique developed by Reegen et al. (2006) 
for MOST Fabry Imaging photometry. To define a fixed 
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Figure 1. The light curve of AQ Leo obtained by MOST. The 
entire time series (upper panel), as well as an expanded view of 
3.7 days of the light curve (lower panel). 

aperture mask by which pixels illuminated by the star and 
by the background can be identified, the PSFs of the star 
on all subraster images must be aligned. This is not the case 
in the raw images due to minor pointing errors during the 
observing run. As a first step, the centroid position of the 
stellar PSF for each image is computed. Subsequently, all 
other pixel intensities relative to the centroid coordinates 
are computed by 2-dimensional linear interpolation. Finally, 
all pixels of each image are shifted by integer values in x 
and y to a pre-defined centre position. Since this shifting 
process decreases the resulting image size common to all ex- 
posures, frames suffering from pointing errors larger than 
3 pixels (about 9 arcsec) from the subraster centre are re- 
jected. (This pointing threshold eliminated 7% of all the ex- 
posures in the raw AQ Leo time series.) Once all the images 
are aligned, the decorrelation can be computed to reduce 
the background effects of modulated stray light. 

Both reductions led to consistent light curves and fre- 
quency identifications. The IfA light curve has higher point- 
to-point scatter than the UBC reduction, and the UBC re- 
duction is presented here as the AQ Leo light curve. How- 
ever, both light curves were searched independently by the 
frequency analysis techniques described in the next section, 
and the results are consistent for both reductions. The final 
duty cycle of the reduced photometry is 82%. The reduced 
light curve is presented in Figure Q] The AQ Leo raw data 
and reduced light curve are available in the MOST Public 
Data Archive, accessed via the Science page of the MOST 
web site: www.astro.ubc.ca/MOST. 



3 FREQUENCY ANALYSIS 

Frequency analyses of the MOST AQ Leo light curve were 
performed independently with three different routines: Pe- 
riod04 (Lenz & Breger 2005), SigSpec (Reegen 2007) and 
nonlinear least squares fitting and bootstrapping (Rowe et 
al. 2006b; see also Cameron et al. 2006) and the results com- 



MOST photometry of the RRd Lyrae variable AQ Leo 3 



pared to ensure that only significant intrinsic stellar frequen- 
cies were identified. 



3.1 Period04 

Period04 (Lenz & Breger 2005) is a routine which applies 
a single-frequency power spectrum and simultaneous multi- 
frequency sine-wave fitting. The program includes advanced 
options such as the calculation of optimal light curve fits 
for multi-periodic signals including harmonics, combination 
frequencies, and combs of equally spaced frequencies. In 
this analysis, we identified frequencies through successive 
prewhitening. As each new significant frequency was added 
to the fit, the optimum values of first overtone frequency fi 
and the fundamental /o, and their amplitudes and phases, 
were determined by minimising the residual rms error of 
the fit. After no additional frequencies were detected in the 
data, we applied a simultaneous least-squares fit of all the 
independent frequencies and their significant linear combi- 
nations, as well as all amplitudes and phases. This included 
a final optimisation of the first overtone and radial funda- 
mental frequencies. 

Besides determining the uncertainties from the error 
matrix of a least-squares calculation and/or from analyti- 
cally derived formulae, Period04 also provides an interface 
for estimating the uncertainties of the fit parameters by 
means of Monte Carlo simulations. A set of time series is 
generated for which the times are the same as for the orig- 
inal time series, and for which the magnitudes (intensities) 
are calculated from those predicted by the last fit plus Gaus- 
sian noise. A least-squares calculation is made for every time 
string, and the resulting uncertainties are based on the dis- 
tribution of the fit parameters. 



3.2 Nonlinear least-squares fitting and 
bootstrapping 

The nonlinear least-squares approach we used for determin- 
ing the sinusoid parameters from the AQ Leo time series 
follows the philosophy of Period04 (Lenz & Breger 2005) 
described above. Discrete Fourier Transforms (DFTs) are 
calculated and an updated fit is subtracted from the data 
successively, until there is no meaningful change in the fit 
residuals. The data are fit using an equation of the form 

mag = A + Aj cos(27r/jt - 6j), (1) 

3=1,™ 

where Ao is a linear offset and fj, Aj and 8j are the fre- 
quency, amplitude and phase for each successive peak found 
in the amplitude spectrum. The uncertainties calculated for 
time series parameters derived from nonlinear least-squares 
fitting depend on the noise of the data (which may be a 
combination of instrumental and random processes) and on 
the time sampling. 

It is also known that fitted phase and frequency param- 
eters are correlated, leading to underestimated uncertain- 
ties in these parameters when calculated from a covariance 
matrix (see, e.g., Montgomery & O'Donoghue 1999). The 
"bootstrap" (see Wall & Jenkins 2003) is a very effective way 
to assess the uncertainties in these fitted parameters. The 



procedure has recently been used in a number of MOST ap- 
plications (see, e.g., Cameron et al. 2006, Rowe et al. 2006b, 
Saio et al. 2006). Clement et al. (1992) also used a boot- 
strap to estimate the uncertainties in Fourier parameters 
they derived for RR Lyrae stars. Bootstrapping produces a 
distribution for each calculated parameter by constructing 
a large number of light curves from the original data. No as- 
sumptions must be made about noise properties of the data 
and individual photometric errors are not required for the 
calculations. Each new light curve is assembled by randomly 
selecting N points from the original light curve (also con- 
taining N points) with the possibility of replacement. The 
new synthetic light curves preserve the noise properties of 
the original data. The fit is repeated for each new light curve, 
eventually building distributions in each of the fit parame- 
ters. We then estimate the la error bars from the analytic 
expression for the standard deviation of each distribution 
under the assumption that they are normally distributed. 
Each distribution is checked to ensure this assumption is 
valid. The uncertainties are listed in Table 1 and give con- 
sistent answers with SigSpec, as described below. We point 
out that bootstrapping only estimates the uncertainties in 
parameters. It does not refine the parameter values or assign 
significances to them. 



3.3 SigSpec and background correction 

A third independent frequency search was performed, where 
the statistical significances of peaks in the DFT of the AQ 
Leo light curve were estimated through SigSpec (Reegen 
2007). This routine calculates an unbiased False Alarm 
Probability (based on white noise) associated with the am- 
plitude of a peak in Fourier space. It produces a significance 
spectrum containing all frequencies above a certain thresh- 
old in significance, as determined by successive prewhiten- 
ing. 

The analyses above identify frequencies present in the 
data with high statistical significance and signal-to-noise, 
but some of these may be due to background variations in 
the data, not intrinsic to the star. The frequencies of some 
background artifacts in MOST data are known; e.g., the or- 
bital frequency of the satellite and its harmonics, at which 
stray light modulation is observed. We carried out tests to 
distinguish periodic effects in the background from intrinsic 
stellar signal. We did this by analysing not only the vari- 
ation of the mean integrated intensity within an aperture 
surrounding the star ("stellar signal", which includes the 
background in those pixels), but also the mean integrated 
intensity outside a larger aperture which includes the out- 
ermost pixels of the CCD subraster ("background signal"). 
Since our three independent frequency analyses yield con- 
sistent results, we chose to use the list of frequencies, ampli- 
tudes and phases from SigSpec to compare the significance 
spectra of the stellar signal light curve and the background 
signal light curve. 

First, peaks from the stellar signal DFT and the back- 
ground DFT are flagged if their frequencies agree within 
the uncertainty. Based on extensive numerical simulations, 
Kallinge, Reegen & Weiss (2007) found the upper limit for 
the frequency uncertainty in a time series to be 
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AT*y/aig(AJ,0) 

where AT is the total data set length in days and sig(A, f, 6) 
is the spectral significance for a given amplitude, frequency, 
and phase, as defined in Reegen (2007). This upper limit 
on frequency uncertainty is more conservative (« 4 times 
larger) than the value for the frequency error estimated by 
Montgomery & O'Donoghue (1999), based on an analytical 
solution for the error of a least-squares fit of a sinusoidal 
signal. If more than one coincident background peak is found 
around a given stellar peak within the frequency resolution, 
then the background peak with the highest significance is 
used for comparison to obtain the most conservative (and 
presumably the safest) solution. 

Next, assuming stray-light induced artifacts are addi- 
tive in terms of their intensities, the peak amplitudes found 
to be significant in the background are rescaled from the 
mean background signal level to the mean stellar signal level 
in order to make the relative amplitudes (and hence the sig- 
nificances) comparable. Now, with the amplitude and spec- 
tral significance of a stellar Fourier peak and the rescaled 
amplitude and spectral significance of its coincident back- 
ground peak, the "conditional probability" that the stellar 
peak is not induced by background is calculated. The False 
Alarm Probability is transformed into a conditional spec- 
tral significance, and only those stellar peaks for which this 
is above 5.46 (corresponding roughly to an amplitude ra- 
tio of 4 for the target peak and the rescaled background 
peak) survive the cut. A detailed publication on this method, 
dubbed CINDERELLA (Comparison of INDEpendent REL- 
ative Least-squares Amplitudes), is in preparation (Reegen 
et al. 2007). 

Uncertainties in frequency, amplitude and phase were 
first estimated through the relations of Montgomery & 
O'Donoghue (1999) and also approximated according to 
Kallinger et al. (2007). The resulting uncertainties are con- 
sistent with the bootstrap values described above, except 
where two frequencies are close to the resolution limit. In 
these cases, which apply to only 4 pairs of identified fre- 
quencies, the least-squares fitting can be perturbed, adding 
further small systematic errors added to the fitted values. 

3.4 The largest intrinsic variations observed by 
MOST - an unexpected challenge 

MOST was designed to detect and characterise very low- 
amplitude stellar oscillations, down to a few /xmag. AQ Leo 
is an exception to the MOST target list, varying over a range 
of a few x 0.1 mag. A common MOST data reduction pro- 
cedure of correcting the data is to apply a running mean to 
reduce the MOST-orbit-modulated stray light signal. In this 
case, however, the stellar variability is large enough that this 
introduces a modulation because the relative contribution of 
the stray light variation changes with the mean brightness 
of the star. By fitting sky values versus instrumental mag- 
nitude, the effects of straylight can be largely removed, but 
the intrinsic stellar variation serves to modulate the relative 
contribution of the stray light to the light curve. This intro- 
duces a small artifact signal sidelobe mirroring the MOST 
orbital frequency but reduced in amplitude compared to the 
intrinsic stellar signal by a factor of about 100. For exam- 



ple, the frequency /i in Table [T] will be reproduced at a 
frequency equal to the MOST orbital frequency minus /i, 
but at a much lower amplitude. 

We tested this through simulations by mimicking the 
raw stellar photometry by adding a constant offset to the 
background light curve and subsequently introducing the 
undeniably intrinsic stellar variability represented by the 
linear combinations of the first overtone /i and fundamen- 
tal mode fo of AQ Leo. As suspected, in addition to the 
background frequencies and the aforementioned linear com- 
binations, peaks due to beat frequencies among f\, fo and 
the MOST orbital period appear in the amplitude spectrum. 
These frequencies are recognised and rejected in the analyses 
of versions of the light curve reduced with a running mean, 
and are not seen in reductions which do not apply this filter- 
ing technique. The MOST team is now alert to this effect for 
any possible future observations of large-amplitude variable 
stars. 



3.5 Intrinsic frequencies in AQ Leo 

As a result of the described procedures, we identified 42 
frequencies in the AQ Leo time series intrinsic to the star. 
All the approaches described above led to this consistent 
solution. The two dominant frequencies are the first radial 
overtone f\ and the fundamental radial mode fo- Of the re- 
maining frequencies, 32 are linear combinations of these two 
up to the 6th order (see the discussion in the next section). 
The uncertainties of the frequencies, amplitudes and phases 
were estimated through extensive Monte Carlo simulations 
as well as through bootstrapping, as described above. The 
resulting uncertainties are comparable with those obtained 
using the formulation of Montgomery & O'Donoghue (1999) . 
All 42 frequencies pass the conditional spectral significance 
test and cannot be reproduced from the background time 
series which contains solely instrumental signal and stray 
light. 

The 42 frequencies include 11 which are not linear com- 
binations of fo and /i. These are 1.96 d _1 , one which cor- 
responds to its first harmonic at 3.92 d _1 (and which has 
a higher amplitude), and their linear combinations with fo 
and /i. We were very cautious about these frequencies be- 
cause of their proximity to integer numbers of cycles per 
day. The nearly continuous MOST time series does not suf- 
fer from cycle/day aliases found in groundbased single-site 
data due to to night/day gaps, but there are subtle modu- 
lations of stray light at 1, 2, 3 and 4 d" 1 due to MOST's 
Sun-synchronous orbit (which brings the satellite to nearly 
the same point over the Earth after 1 day). However, our 
analysis identified and rejected those artifacts at 2.00 and 
4.00 d _1 present in the background measurements and did 
not reject the frequencies at 1.96 and 3.92 d , which differ 
from potential artifacts by more than the frequency resolu- 
tion of the data set. Indeed, these frequencies are retained 
with high significance and signal-to-noise. To stress their im- 
portance among the set of 42 frequencies, we hereafter refer 
to them as fi = 1.96 d _1 and fu = 3.92 d" 1 . 

The frequencies, their amplitudes, phases and uncer- 
tainties are listed in Table[T] The amplitude spectrum show- 
ing these frequencies is plotted in Figure [2] 
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Table 1. Frequencies identified in AQ Leo from the MOST photometry by SigSpec, which are consistent with the other frequency 
analysis techniques described in the text. Amplitudes, phases, significances, and uncertainties (derived by bootstrapping) arc also listed. 
The phases correspond to epoch JD2000 = JD2451545. 

id: identification, /; frequency [d" 1 ], erf. error in frequency [d — 1 ] (lcr), A: amplitude [mmag], a a'- error in amplitude [mmag] (lcr), sig.: 
spectral significance, 8: phase, erg: error in phase (lcr), comb. err.: deviation / — f c from the numerical value for combinations 

fc = afo + bfx + cfi 
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4 DISCUSSION 

4.1 Combination frequencies of the two dominant 
modes 

The MOST photometry has revealed an unprecedented 
number of frequencies (42) in AQ Leo. Jerzykiewicz & Wen- 
zel (1977) (hereafter JW77) first detected the first radial 
overtone and fundamental radial modes in this star. They 
modeled their light curve with a double-harmonic series and 
by fitting different numbers of higher-order terms until the 
standard deviation of the fitted curve was near the lower 
limit of the mean error of a single observation. JW77 arrived 



at a fit with 14 frequencies, of which 12 are combination fre- 
quencies of the two dominant modes. Of our 42 frequencies, 
32 are linear combinations of those same modes. 

The numerical combinations can be obtained arithmeti- 
cally in Table[T]but are much easier to see when the frequen- 
cies are plotted in " echelle" diagrams like Figure [3] These 
diagrams have the same form as the echelle diagrams used 
in asteroseismology to highlight the asymptotic distribution 
of nonradial p-modes. The frequencies are folded modulo 
a common frequency spacing on the x-axis and plotted at 
their original values on the y-axis. In the case of the MOST 
AQ Leo frequencies (and for RRd frequencies in general), 
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Figure 2. Upper panel : The amplitude spectrum of the AQ Leo 
light curve. Most of the linear combinations of fo and /i can 
be identified by eye. Lower Panel : The amplitude spectrum for 
frequencies identified by SigSpec. The inset is a close-up of the 
low frequency range containing fi and fa. 



the frequencies align clearly in such a diagram due to their 
arithmetic relationships. 

In Figure[3j the points on the right-hand side of the dia- 
gram and climbing upwards to the left with regular "jumps" 
are the linear combinations of fo and /i, or i-fi + j-fo- 
(For this diagram, they are written in the practical form 
of k-fi + l-fo + m(fi — fo), where i,j and k,l,m are inte- 
ger numbers.) The "jumps" in the diagram occur when the 
linear combinations go up by one order, and are due to the 
simple fact that /imod(fi — fo) = fomod(fi — fo) 7^ 0. Hence, 
for every additional fo or f\ in a linear combination (besides 
multiples of (/1 — fo)), the "folded" (modulo) value changes. 
The change corresponds to 8f = 3-/i - 4-/ = 0.03879 d" 1 . 
The points at almost the same y-value in the diagram show 
that, from this MOST data set, we were able to separate dif- 
ferent linear combinations of fo and /1 as close as 0.03879 
d" 1 to one another. The frequencies which fall in the left 
bottom corner of Figure[3]are associated with the additional 
frequencies at 1.96 and 3.92 d _1 . Note that those frequencies 
also follow a pattern of regular "jumps", since we have de- 
tected linear combinations of these frequencies with fo and 
/1. We stress that the regular patterns are merely a conse- 
quence of the arithmetic relations between the frequencies 
in the linear combinations. 
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Figure 3. Graphical presentations of the linear combinations of 
the two dominant modes in AQ Leo, in the form of "echelle" 
diagrams. The observed frequencies folded at the beat frequency 
/1 — /o = 0.61925 d" 1 show differences of Sf = 3 /i — 4-/o in a 
regular pattern alternating by fo and f±. 



The MOST photometry of AQ Leo opens the possibility 
of a detailed analysis of the properties of the combination 
frequencies resulting from nonlinear coupling, which can be 
used to test hydrodynamical models of RRd stars (see Sec- 
tion 4.2). 

The phase diagrams of the two dominant modes are 
shown in Figure [4] in which the linear combinations and 
frequencies unrelated to the two dominant modes we find in 
the fit (Table [T]) have been subtracted. In AQ Leo, the first 
overtone mode is the dominant one. It is characterised by 
stronger nonlinear behaviour which can be seen by the more 
pronounced nonsinusoidal shape of its light curve in Figure 



SI 



4.2 Amplitude ratios and comparison to theory 

Early in the study of RRd stars, it was believed that they 
were in the process of switching from the fundamental ra- 
dial mode to the first overtone or vice versa. Cox et al. 
(1980) concluded that the double- mode behaviour of AQ Leo 
could be explained by mode switching during its blueward or 
redward horizontal branch evolution. This mode-switching 
phase is very brief compared to the star's duration as an RR 
Lyrae variable - lasting only a few thousand years (Cox et 
al. 1980, 1983). This is too short to account for the observed 
numbers of RRd stars in clusters like M15 (see Paparo et al. 
1998), which suggests that the RRd stage must on average 
last longer than the theoretically predicted mode-switching 
time scales. Feuchtinger (1998) was the first to find stable 
(i.e., not resulting from mode switching) double-mode be- 
haviour in RR Lyrae models by including a time-dependent 
mixing length model for convection. Subsequently, convec- 
tion was established by Kollath et al. (2002) as the key player 
in double-mode pulsation in Cepheid and RR Lyrae models 

Even so, there is some observational evidence that the 
transition from a monoperiodic RR Lyrae star (RRc: red- 
ward, or RRab: blueward) to a double-mode RR Lyrae 
star can take place over a very short timescale. Clement & 
Goranskij (1999) observed a rapid mode change from RRab 
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Figure 4. The deconvolved residuals of the MOST light curve 
of AQ Leo (after removal of combination frequencies and the ad- 
ditional frequencies) folded with the fundamental radial mode 
period Po (upper panel) and with the first overtone radial mode 
period Pi (lower panel). 



to RRd in the star V79 in the cluster M3 in less than a 
year. Buchler & Kollath (2002) have shown that bifurcation 
points (characterized by changes in pulsation behaviour) ex- 
ist where pulsation is no longer independent from stellar 
evolution. Taking this into account, the number of observed 
double-mode pulsators suggests that these oscillations are 
stable and that rapid mode switching is a relatively rare 
effect. Still, the true nature of double-mode pulsation in in- 
dividual stars remains elusive from the observers' point of 
view. 

If the kinetic energy switching times are only a few 
hundred years or less and the relative amplitude change 
timescale less than a thousand years, the relative amplitudes 
of the two modes might be seen to change over timescales of 
decades. The amplitude ratio of the first overtone to funda- 
mental mode in AQ Leo in 1973-1975 measured by JW77 was 
Ai/Aq ~ 2. The amplitude ratio measured by MOST in 2005 
is Ai/Ao — 1.8. Unfortunately, it is risky to compare our ob- 
served amplitudes and phases to those tabulated by JW77, 
because the MOST photometry was obtained through a wide 
nonstandard passband while the JW77 data were obtained 
through a Johnson B filter. 

A qualitative comparison to a RRd model (see Kovacs 
& Buchler 1993) is presented in Figures [S] and [S] The strik- 
ing similarities to the model results in the MOST light curve 
and resulting amplitude spectrum shows that current space- 
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Figure 5. Comparing light curves from a computer and from 
space. A segment of the theoretical light curve of a purely ra- 
diative RRd star, taken from Kovacs & Buchler (1993) covering 
10 cycles of the fundamental mode period in their model (upper 
panel - reproduced by permission of the AAS). A segment of the 
MOST light curve of AQ Leo spanning the same amount of time, 
shown at the same scale (lower panel). It is now possible to ob- 
tain data which is sampled as thoroughly and with noise levels 
comparable to time-resolved model calculations. 



based observations are able to reach the quality of numerical 
simulations. Admittedly, the frequencies of the fundamen- 
tal and the first overtone themselves don't match and we 
do not propose this model to be an accurate description of 
AQ Leo. The model was calculated under the unrealistic as- 
sumption of a purely radiative stellar envelope. By today's 
standards, more sophisticated treatments of the hydrody- 
namic equations (including convection) in the calculations 
are more reasonable (Kollath & Buchler 2001). Still, the pre- 
dicted relative amplitudes of the various linear combinations 
of fo and /i behave in close agreement with our observations 
up to very high order. As such, observations of the quality 
presented in this paper should be considered to be a tool of 
validation for the non-linear pulsational behaviour of mod- 
ern RRd models. The observed frequencies, their amplitudes 
and phases can even be a testing ground for the development 
of hydrodynamical models of RRd stars including nonradial 
pulsations. 



4.3 Period changes? 

Basic pulsation physics predicts that an increase in the fun- 
damental mode period would be accompanied by an increase 
in the first overtone mode period, reflecting an increase in 
the mean stellar density. Paparo et al. (1998) studied period 
changes in RRd stars in M15, and found that the rate of pe- 
riod change for the fundamental radial mode is significantly 
larger than for the first overtone. For most of the stars in 
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Figure 6. Amplitude spectra of the model light curve from Fig. [5] 
taken from Kovacs & Buchler (1993) (upper panel - reproduced by 
permission of the AAS) and the MOST photometry in the same 
frequency range (lower panel). The relative amplitudes compare 
quite well to our findings. 

their sample, the ratio Po/Pi is increasing. They found the 
period changes of the fundamental modes of the RRd stars 
were mostly negative, while for the first overtone, both pos- 
itive and negative changes were seen. The most striking re- 
sult by Paparo et al. (1998) is the measurement of period 
changes of different sign for the fundamental and first over- 
tone modes. This suggests, as for the 8 Scuti pulsators, that 
the period changes may not be simply linked to an RR Lyrae 
star's evolution in the HR diagram. 

The uncertainties in the frequencies in Table [T] are 
slightly too large, given that the MOST time series spans 
only about 34 days, to measure reliably the predicted evolu- 
tionary period changes in the 30 years between the MOST 
and JW77 observing runs. But Jerzykiewicz, Schult & Wen- 
zel (1982) (hereafter JSW82) deduced that there was a sud- 
den increase of the period of the first overtone of AQ Leo in 
the early 1970s, by comparing re-evaluated JW77 photome- 
try to archived data dating back to 1929. It is possible that 
there have been other period changes in AQ Leo unrelated 
to the evolution of the mean density of the star. 

Therefore, we compared the periods from the MOST 
photometry to the values reported in JSW82. (We cannot 
perform a meaningful O — C analysis due to the complexities 
of the light curves.) Our results yield evidence for frequency 
changes in AQ Leo since the JW77 observations outside the 
la uncertainties: fundamental mode, 

A/ = -0.00005 ± 0.00002 d" 1 

and first overtone, 

A/i = +0.00007 ± 0.00001 eT 1 . 

While the change of the fundamental frequency is consis- 
tent with zero within the 3a confidence limits, the non-zero 
positive change of the first overtone is highly significant. 
Interestingly, the value of /i we obtained from the MOST 
data matches that proposed in JSW82 to satisfy an O — C 
diagram covering the data obtained from 1929 up to right 



before the abrupt period change in the 1970s. Our results fit 
their findings to within the la uncertainties. This hints that 
the period change in the 1970's was a short-lived event. 

Consequently, notwithstanding the higher uncertain- 
ties, we obtain a slightly different period ratio than that 
of JSW82, namely P 1 /P = 0.74602 ± 0.00002 compared 
to 0.746063 ± 0.000003. This does not significantly disagree 
with the theoretical result of Cox et al. (1980), who cal- 
culated that the fundamental mode period and the period 
ratio of 0.746 can be explained if AQ Leo has a homoge- 
neous composition typical of Population II and a mass of 
0.65 Mq. More recent calculations (Kollath, private commu- 
nication) indicate, however, that a mass > 0.7 Mq is likely to 
be more accurate . When compared to the results of Szabo, 
Kollath & Buchler (2004), our values match the sequence 
with M — 0.71 Mq, L = 60 Lq and very low metallicity. 



4.4 New modes in AQ Leo? 

RR Lyrae stars have long been considered prototypes of radi- 
ally pulsating stars. Two types of RR Lyrae stars show mul- 
tiperiodic behaviour: Blazhko stars and double-mode stars, 
like AQ Leo. The former are characterised by the Blazhko 
Effect (see Kolenberg 2004 and Kovacs 2001), which is a 
periodic modulation of the amplitude and/or phase of the 
light curve on timescales of a few up to hundreds of days. 
In the Fourier spectra of such stars, additional frequencies 
occur close to the dominant frequency which is identified 
with a radial mode. If these frequencies correspond to pul- 
sation modes in the star, they must be nonradial modes. Van 
Hoolst, Dziembowski & Kawaler (1998) showed that there is 
a dense spectrum of nonradial modes which have lower mo- 
ments of inertia (and hence are more easily excited) in the 
vicinity of the radial mode frequencies in RR Lyrae stars. 

The MOST observations of AQ Leo are the first in 
which additional frequencies besides the fundamental ra- 
dial mode, the first overtone radial mode and their linear 
combinations are detected. The two additional frequencies 
are fi ~ 1.96 dT 1 and fu ~ 3.92 d~ . (Linear combinations 
of these frequencies with the known radial modes are also 
found in the data, see also Figure [3]) The amplitude ratio 
of fi to fu is 0.64 ± 0.08. 

Judging from the period ratios alone, fu (~ Ifi) could 
be the third radial overtone (Zoltan Kollath, private com- 
munication). If it is a strange mode (essentially an acoustic 
surface mode; see Buchler, Yecko & Kollath 1997), /, could 
then be a nonradial mode in a resonant (2:1) interaction with 
the third overtone. This scenario has been considered theo- 
retically in studies of pulsational resonance (Buchler, Goupil 
& Hansen 1997), supporting our argument. If we are indeed 
witnessing resonance of a non-radial and a radial mode, it is 
the higher frequency which must correspond to t = due to 
parity considerations. The frequencies match the resonance 
within the uncertainties (see Table 1). In case of a mismatch, 
the resonance cannot be locked. The fact that we see sig- 
nificant amplitudes at both frequencies may suggest that 
AQ Leo is in a fleeting state of mode transition. Previous 
calculations by Buchler, Yecko & Kollath (1997) predicted 
that if strange modes are present, they should correspond 
to higher overtones than the third. Our new frequencies are 
more consistent with the third overtone, so models matched 
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more closely to AQ Leo are necessary to compare with the 
MOST observations. 

The frequency fi - close to the radial fundamental mode 
- is also reminiscent of the additional frequencies seen in 
Blazhko stars. While the nature of the additional frequencies 
in AQ Leo, as well as in Blazhko stars (see, e.g., Kolenberg 
et al. 2003), must still be thoroughly investigated, there is 
increasing evidence that some RR Lyrae stars pulsate "be- 
yond radial modes" (e.g., Clement & Rowe 2000; Kiss et 
al. 1999). Our observations represent the first detection of 
additional modes in an RRd star. 



5 WHERE DO WE GO FROM HERE? 

Nearly continuous MOST photometry of AQ Leo spanning 
about 34.4 days and with high photometric precision have 
provided (1) accurate Fourier components of the known ra- 
dial fundamental and first overtone modes in this star; (2) 
an unprecedented number of 32 combination frequencies of 
these modes; (3) evidence for non-evolutionary changes in 
the periods of the fundamental and (even more strongly) 
the first overtone mode; and (4) additional frequencies which 
can only be plausibly explained by other modes which have 
never been seen before in AQ Leo, or in any RRd star until 
now. What are the next steps forward? 

Observationally, it will be very difficult - perhaps im- 
possible - to obtain a light curve of longer time coverage 
and better precision than the MOST photometry presented 
here in the near future. The star will not be accessible to 
the CoRoT and Kepler space missions, nor will groundbased 
photometric campaigns reach what MOST has achieved for 
AQ Leo for some time to come. If AQ Leo is reobserved 
by MOST in the coming years, it could provide a reliable 
measurement of the change in amplitude ratio of the radial 
fundamental and overtone modes in the MOST bandpass, to 
test theories of mode switching and evolution in this star. 
Also, further evidence for the period changes may be estab- 
lished. 

New MOST photometry may lead to the detection of 
even more frequencies (depending on the behaviour of AQ 
Leo) but it will not be possible to improve the frequency res- 
olution of the current data set dramatically. Time-resolved 
high-resolution spectroscopy to look for evidence of nonra- 
dial pulsation in the spectral line profile shapes would be 
highly desirable for this star. Unfortunately, AQ Leo is quite 
faint (V ~ 12.6), so spectroscopy of sufficient wavelength 
resolution and high signal-to-noise are beyond reach of in- 
struments which could monitor the star at a sufficient ca- 
dence for many weeks. 

The immediate path to better understanding of AQ Leo 
is theoretical modeling of the frequency spectrum identified 
in the existing MOST photometry. The potential of the com- 
bination frequencies to test hydrodynamical models can now 
be fully exploited thanks to the very precise MOST data. 
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